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We have analyzed the symmetry properties and the ground state of an orbital Hubbard model 
with two orbital flavors, describing a partly filled spin-polarized e g band on a cubic lattice, as in 
ferromagnetic manganites. We demonstrate that the off-diagonal hopping responsible for transitions 
between x 2 — y 2 and Sz 2 — r 2 orbitals, and the absence of SU(2) invariance in orbital space, have 
important implications. One finds that superexchange contributes in all orbital ordered states, the 
Nagaoka theorem does not apply, and the kinetic energy is much enhanced as compared with the spin 
case. Therefore, orbital ordered states are harder to stabilize in the Hartree-Fock approximation 
(HFA), and the onset of a uniform ferro-orbital polarization and antiferro-orbital instability are 
similar to each other, unlike in spin case. Next we formulate a cubic (gauge) invariant slave boson 
approach using the orbitals with complex coefficients. In the mean-field approximation it leads 
to the renormalization of the kinetic energy, and provides a reliable estimate for the ground state 
energy of the disordered state. Using this approach one finds that the HFA fails qualitatively in the 
regime of large Coulomb repulsion U — ► oo — the orbital order is unstable, and instead a strongly 
correlated orbital liquid with disordered orbitals is realized at any electron filling. 
[Published in: Phys. Rev. B 71, 144422 (2005).] 
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I. INTRODUCTION 



In recent years there has been renewed interest in 
orbital degrees of freedom in Mott insulators** Typi- 
cally, Mott insulators are stoichiometric, i.e., undoped, 
oxides (or sulfides) in which the strong on-site interor- 
bital Coulomb repulsion U on the transition-metal ions 
dominates over the kinetic energy driven by the electron 
hopping i, and eliminates charge fluctuations. At en- 
ergies well below U one is then left with effective low- 
energy interactions cx t 2 /U of superexchange (SE) type. 
In many cases these are purely magnetic interactions be- 
tween the spins on the metal ions, leading to the familiar 
spin models, such as the Heisenberg model. However, 
when the electrons occupy partly-filled degenerate e g or 
t2 g orbitals, such as in the perovskites KCUF3, LaMnC>3, 
LaTi03, and LaVC-3, the orbital degrees of freedom be- 
come equally important as the spin ones, and it is there- 
fore necessary to treat both of them on equal footing. In 
such cases the SE is described by so-called spin-orbital 
models?^ and the SE interactions are typically strongly 
frustrated, even on a cubic latticed In spin-orbital models 
the quantum effects are particularly strong — the quan- 
tum fluctuations are enhanced, and might even desta- 
bilize the long-range magnetic order, leading to a spin 
liquid state, possibly realized in LiNi02^ The opposite 
situation, that an (isotropic or anisotropic) orbital liq- 
uid (OL) is stabilized and coexists with long-range spin 
order, was pointed out recently for t2 g Mott-Hubbard 



systems^ By contrast, in undoped e g systems, such as 
KCuF 3 (Ref. 0) and LaMn0 3 (Ref. %, the quantum 
phenomena are partly quenched and the SE favors alter- 
nating orbital (AO) order which coexists with antiferro- 
magnetic (AF) spin order. 

An issue of considerable interest is how such systems, 
characterized by the presence of orbital degrees of free- 
dom, behave under doping, and in particular how this 
compares with the more familiar behavior of doped spin 
systems. In this paper we address this issue by study- 
ing a generic model of correlated e g electrons with two 
orbital flavors, described by a pseudospin T — 1/2 in 
the orbital Hilbert space, and consider its relation to the 
standard (spin) Hubbard model for electrons with spin 
S = 1/2. So we introduce the e g -orbital Hubbard model 
and investigate: (i) in what respect long-range order in 
such an orbital system is different from that in the analo- 
gous spin system, and (ii) whether the orbitals may order 
when U is large, or rather form a disordered OL. These 
questions are of fundamental nature and our main aim in 
addressing them is to uncover and elucidate the physical 
mechanisms which operate in the e g band and are typical 
for orbital degeneracy, in particular by contrasting them 
with those known to operate in spin systems. 

The present problem is closely related to the physi- 
cal properties of the colossal magnetoresistance (CMR) 
manganites,- where the well-known mechanism of dou- 
ble exchange introduced by Zcncr 8 is responsible for the 
metallic ferromagnetic (FM) phase at finite doping, in 
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which the spins of the e g electrons are fully polarized. 
The model that we will investigate here covers only the 
case of the FM phase, thus neglecting the competition 
of the double-exchange mechanism with the spin AF SE, 
and the resulting dependence of the hopping amplitude 
on the actual spin states at two neighboring sites£ Even 
when one limits oneself to the FM phase, a sequence 
of orbital-ordered phases may be expected^ and each 
of them would break cubic symmetry, contrary to what 
is observed in the magnetic properties of the metallic 
FM phase. The analysis of the present paper, mak- 
ing an extensive use of the auxiliary particle method in 
the strongly correlated regime^ provides the basis for a 
proper treatment of this problem-^ which enables one to 
understand the persistence of cubic symmetry and more 
in particular why the magnon stiffness constant increases 
with hole doping.— 

This paper is organized as follows. In Section [IT] we 
introduce the orbital Hubbard model for spin-polarized 
e g electrons at orbital degeneracy, and discuss its sym- 
metry properties. We show that the cubic symmetry of 
the hopping may be better appreciated when a particular 
basis consisting of two orbitals with complex coefficients 
is used. Next we analyze in Sec. IIIII the possible orbital 
ordered phases at U — oo and compare their densities of 
states and total energies derived within the slave fcrmion 
formalism. Such phases follow from the instabilities to- 
wards orbital-ordered states obtained within the Hartree- 
Fock (HF) approximation (Sec. IIV|) . and we show that 
such instabilities and the properties of the ordered phases 
at finite U, related to the SE, are here quite different from 
those known from the spin Hubbard model. In Sec. [V] 
we introduce the cubic invariant slave boson approach 
and use the mean-field approximation to analyze the dis- 
ordered orbital liquid state. Within a generalization of 
the Kotliar-Ruckenstein^ (KR) approach to the present 
orbital problem, we give reasons why the orbital ordered 
states are unstable against the OL disordered state when 
one goes beyond the HF approximation. The paper is 
concluded in Sec. IVII by pointing out the implications of 
our results for the physical properties of the CMR man- 
ganites. 



II. ORBITAL HUBBARD MODEL 

A. The Hamiltonian and its symmetry properties 

We consider spinless e g electrons on a cubic lattice with 
kinetic energy 

a (ij)\\a 

where hopping with amplitude —t between sites i and j 
occurs only for a pair of directional orbitals |£ a ) oriented 
along the bond (ij) direction, i.e., |( Q ) oc 3x 2 — r 2 , 3y 2 — 
r 2 , and 3z 2 — r 2 , when the bond (ij) is along the cubic axis 



a = a, 6, and c, respectively. We will similarly denote by 
\£ a ) the orbital which is orthogonal to |£ a ) and is oriented 
perpendicular to the bond (ij), i.e., |£ Q ) oc y 2 —z 2 , z 2 ~x 2 , 
and x 2 — y 2 , for a bond (ij) along the axis a — a, b, and c, 
respectively. While such a choice of basis, that depends 
on the bond direction under consideration, is convenient 
for writing down the kinetic energy, one cannot avoid to 
choose a particular orthogonal basis for the two orbital 
flavors as soon as one wants to introduce a Hubbard term 
to describe the local electron interactions. The usual 
choice is to take 

|z)^(3z 2 -r 2 ), | x )^^(x 2 -y 2 ), (2.2) 

called real orbitals. However, because this basis is the 
natural one only for the bonds parallel to the c axis but 
not for those in the (a, b) plane, the kinetic energy then 
takes the forrr*i£i±£ 

H t = [ 3c lx c jx + 4z c jz T V3(cl x c JZ + c\ z c JX )} 

"* E 4 C ;*> ( 2 - 3 ) 

and although this expression is of course cubic invariant, 
the representation (3) of the hopping does not exhibit 
this symmetry but takes a very different appearance de- 
pending on the bond direction. 

We thus prefer to use instead the basis of complex or- 
bitals at each siteii 

|+> = ^(|*>-i|*», h) = ^(k)+^)), (2.4) 

corresponding to "up" and "down" pseudospin flavors, 
with the local pseudospin operators defined as 

-^i = C i+ C i-J = C i- C i+' 

n = h(4+Ci+-4-Ci-) = U n i+- n i-)- ( 2 - 5 ) 

For later reference it is convenient to introduce also elec- 
tron creation operators cf(V>i,#i) which create e g elec- 
trons in orbital coherent states, defined as 

10,) = er^/ 2 cos (*) \i+) + e^ 2 sin (|) |z-)<2.6) 

in analogy with the well-known spin coherent statesJ^ 
The expectation value of the local pseudospin operator in 
the coherent orbital (|2.6I) behaves like a classical vector^ 

(fiilTilfij) = ±(sin7/; l cos6l l ,sinV> l sin6l l ,cosV> l ), ( 2 -7) 

traversing a sphere, with the "equatorial plane" (ipi = 
7r/2) corresponding to real orbitals |Oj(7r/2, 9i)) = \iOi) = 
cos(0j/2)|zz) —sm(8i/2)\ix), and the "poles" (ipi = and 
ipi = 7r) to the complex orbitals \i+) and \i— ). The three 
directional orbitals \iC a ) at site i, associated with the 
three cubic axes (a = a, 6, c), are the real orbitals with 
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9i being equal to d a — — 47r/3, fib — +47r/3, and d c = 0, 
respectively, i.e. 

|<C«) = ^[ e -«-/3|i+> + e+««/ a |i->] 

= cos(i? /2)|ij»> -sin(i? Q /2)|ia;>, (2.8) 

and thus correspond to the pseudospin lying in the equa- 
torial plane and pointing in one of the three equilateral 
"cubic" directions defined by the angles "3 a . 

In the complex-orbital representation (|2.4|) the orbital 
Hubbard model for e g electrons takes the form 

" = -5*E £ [(4 + c j+ +cl Cj _) 

a (ij)\\a 

+l(e- i ^4 + c j _+e +i ^4_ Cj+ j 
+Uj2n l+ n^ 7 (2.9) 

i 

with Xa = +27r/3, \b = — 27r/3, and \ c = 0, and 
where the newly introduced parameter 7, explained be- 
low, takes the value 7=1. The appearance of the phase 
factors e ±lXa is characteristic of the orbital problem — 
they occur because the orbitals have an actual shape in 
real space so that each hopping process depends on the 
bond direction. The form of the interorbital Coulomb 
interaction cx U is invariant under any local basis trans- 
formation to a pair of orthogonal orbitals; it gives an en- 
ergy U either when two real orbitals are simultaneously 
occupied, U Ui X ni Zl or when two complex orbitals are 
occupied, as in Eq. (|2.9|) . 

The representation (|2.9[) has several advantages: (i) It 
displays manifestly the cubic symmetry, since the trans- 
formation Xa — ► Xa + 27r/3 (which amounts in Eq. (|2.9p 
to the cyclic permutation a^6->c-^aof the cubic 
axes) in conjunction with the corresponding phase shift of 
the electron operators cL_ — > c\ ± e ±2lv ^ 3 (which permutes 
the | Ca) -orbitals according to 3x 2 — r 2 — > 3y 2 — r 2 — > 
3z 2 — r 2 — > 3x 2 — r 2 ) leaves the Hamiltonian (|2.9[) in- 
variant, (ii) It exhibits clearly the difference between the 
spin case and the orbital case. In the orbital case there 
is both pseudospin-conserving hopping [the first line in 
Eq. (|2.9p ] and non-pseudospin- conserving hopping [the 
second line in Eq. (|2.9[) ]. whereas in the corresponding 
spin case, i.e., in the standard Hubbard model, there is of 
course only spin-conserving hopping and the second term 
is absent. Thus the present complex-orbital representa- 
tion allows us to introduce the parameter 7 by which 
one can turn the e s -band orbital Hubbard model (7 = 1) 
into what is formally a spin Hubbard model with the 
same hopping amplitudes (7 = 0), interpreting "+" and 
"— " as "spin up" and "spin down". This device makes 
it very easy to recognize the differences in physical be- 
havior between the orbital case and the spin case: the 
parameter 7 will of course show up in all analytical ex- 
pressions below, and one can compare at a glance the re- 
sult for the orbital case (7 = 1) with that for the spin case 
(7 = 0). At the present stage one can already observe 



from Eq. (|2.9[) that there is more kinetic energy available 
per electron in the orbital case, because additional hop- 
ping channels are present. We will see below that this has 
important consequences for the relative stability of vari- 
ous states, (iii) Finally, it shows explicitly that rotational 
SU(2) symmetry for the pseudospins is absent^ 2 - which in 
the complex-orbital representation is immediately obvi- 
ous from the presence of the non-pseudospin-conserving 
hopping term cx 7 in (|2.9|) . Thus the components of 
the total pseudospin operator T — J^Tj, are conserved 
only at 7 = (i.e., [T,H] = 0), while the terms oc 7 in Ti 
commute instead with the staggered pseudospin operator 
Tq = exp(iQ • RJT?, where Q = (tt, tt, tt). 



B. New features compared with the spin case 

It is instructive to follow the changes of the electronic 
structure of the uncorrelated band [i.e. with U = in 
Eq. (|!H))) ] with increasing 7 (0 < 7 < 1). When the 
hopping is only diagonal between pairs of |+) and |— ) 
states at 7 = 0, the pseudospin bands are degenerate, and 
the density of states has the familiar shape obtained for 
a simple cubic lattice, with bandwidth 6t, corresponding 
to the hopping elements of \t [Fig. 0Ja)]. Increasing 7 
removes the degeneracy of the electron bands, and gives 
increasing spectral weight near the band edges without 
modifying the bandwidth. For genuine e g electrons (i.e. 
at 7 = 1) the density of states does not start from zero at 
u> = ±3t, as usually for three-dimensional (3D) lattices, 
but is finite there and has a value close to its average 
over the entire band [Fig. [2d)]. Not only is the spectral 
weight transferred to lower energies, but even the Fermi 
energy at fixed electron density n < 1 decreases with 
increasing 7, as shown on the example of n = 0.7 in Fig. 
[T] Therefore, for a given electron density, at U = the 
kinetic energy of e g electrons (i.e. at 7 = 1) is lower than 
in the corresponding spin case (at 7 = 0). 

Finally, some remarks on the physical interpretation of 
the orbital Hubbard model are in place here. The first of 
them concerns electron spin. As said, the electrons in the 
model are spinless [cf. Eq. (|2.1[) ]. which at first sight may 
seem unphysical. However, such a model is entirely ap- 
propriate for real, i.e. spincarrying e g electrons in a FM 
state, where the spins are fully polarized. This situation 
can be realized in a strong magnetic field, or, as in man- 
ganites, when the double exchange polarizes the <2g core 
spins which in turn polarize the e g band by strong Hund's 
rule coupling. Then the spin degrees of freedom are com- 
pletely frozen out and only the orbital degrees of freedom 
remain and can contribute to the kinetic energy. Actu- 
ally, Eq. (|2.I[) [but with the additional constraint of no 
double occupancy] is precisely the expression for the ki- 
netic energy of the e g band in the metallic ferromagnetic 
phase of the doped manganites Lai-^A^MnOa (with A 
= Sr, Ca,..., and x ~ 0.3) when these are described by an 
extended (orbital-degenerate and large spin) t-J models 
So the U — > 00 limit of the present orbital Hubbard model 



4 



0.6 
0.4 
0.2 
0.0 
0.4 
0.2 
0.0 
0.4 
0.2 
0.0 
0.4 
0.2 
0.0 





-j 1 1 -J 1 1 1 ] 

! \ (a) : 

i \ 


■ 1 1 1 


> i ■ i ■ i ■ 

/\. (b): 

'i l i I i I i 


.■III 




'I'M 


■ I ■ I ■ I ■ 

(d): 



-3-2-10 1 2 3 
co/t 



FIG. 1: Evolution of the density of states N(u>) (in units of 
t — 1) obtained for the tight-binding model (|2.9p at U = 
with increasing off-diagonal hopping jt: (a) 7 = 0, (b) 7 = 
0.5, (c) 7 — l/\/2, (d) 7 = 1; 7 = and 7=1 corresponds 
to the spin and to the orbital Hubbard model, respectively. 
Dashed lines show, for the same electron filling n = 0.7 in all 
cases, the Fermi energy, which decreases with increasing 7 for 
n < 1. 



7 = — 1 to hopping only". Although such £- hopping 
occurs, for instance, in transition metals as a (ddS) el- 
ement, and is symmetry-allowed in the perovskites, it 
cannot occur by the familiar mechanism of two-step hop- 
ping (neither cr-type nor 7r-type) via a 2p orbital on the 
oxygen ion in between two transition metal ions. It is 
therefore generally accepted that in physically relevant 
cases this hopping process is smaller by at least two or- 
ders of magnitude than that between orbitals, and 
thus, to our knowledge, all work on the manganites has 
actually been done assuming pure £- hopping, i.e. 7=1. 
Nevertheless, we will occasionally let 7 vary between 
and 1, not with the intention to suggest that a signif- 
icant strength of ^-hopping is actually physically rele- 
vant, but rather with the purpose of demonstrating how 
the non-pseudospin-conserving hopping affects the phys- 
ical properties of strongly correlated electrons in a partly 
filled band. 

The third remark concerns the difference between real 
and complex orbitals. It is noteworthy that, unlike in 
the spin case, already for an individual site there is no 
spherical symmetry in pseudospin space even at the clas- 
sical level: the directions available to the pseudospin- 
vector are not all physically equivalent. In particular, 
the real orbitals are spatially anisotropic and have a 
nonzero diagonal electric quadrupole moment (EQM), 
(T?) 2 + (Tf) 2 ^ 0, whereas the complex orbitals have 
a cubic shape, with only (Tf) 7^ 0. This difference is of 
course the origin for the hopping Hamiltonian not hav- 
ing SU(2) symmetry. Moreover, as pointed out by Van 
den Brink and Khomskii,-- in a real compound like a 
perovskite the EQM couples directly to the lattice, and 
occupancy of a real orbital would induce a local Jahn- 
Teller (JT) distortion whereas occupancy of a complex 
orbital would noti^i 



III. ORBITAL ORDERED STATES 



A. Uniform and alternating orbital order 



is directly relevant for the physics of the mangan- 
ites, and for this reason we will pay extra attention to 
this limit. 

The second remark concerns the parameter 7. As one 
can readily verify, the kinetic terms in Eq. (12. 9|) with arbi- 
trary 7 are equivalent to the kinetic energy Hamiltonian 



1 

—t , 
2 ^ 



E [( 1 +7)4c c k+( 1 -7)c: 



(2.10) 



which reduces to H t [Eq. (|2.1[) ] for 7=1. So, although 
we have introduced the parameter 7 purely as a formal 
device, it actually describes the relative strength of hop- 
ping between the orbitals perpendicular to a bond, 
and one sees that 7=1 corresponds to "C-hopping only" , 
7 = to "^-hopping and ^-hopping equally strong" 
(equivalent to the spin case as discussed above), and 



Because the electrons interact by the local Coulomb 
interaction U, they are prone to instabilities towards or- 
bital order, similar to the magnetic instabilities towards 
spin order in the spin case^ to which we will compare 
them. At half-filling (n = 1) the simplest possibility to 
reduce the interaction energy cx U would be to polarize 
the system completely into ferro orbital (FO) states, 



I$fo) =Y[c]^,6)\0), 



(3.1) 



with the pseudospin pointing in the same direction at all 
sites. As in the spin case, another possibility is alternat- 
ing orbital (AO) order, 

I*ao) = l[cl(ip A ,8 A ) 11 ctty B ,0 B )|O>, (3.2) 

i£A jeB 
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with orbitals alternating between two sublattices A and 
B which cover a cubic lattice. Depending on whether or- 
bitals alternate in every direction, or whether there are 
lines or planes of ferro orbital order, these states are clas- 
sified as G-type (for spin called Neel states), C-type, or 
A-type AO states. Doubly occupied sites are explicitly 
avoided in all these states. If the band is partly filled 
(n < 1), these ordered states must of course be modified 
to involve a coherent mixture of orbital-polarized occu- 
pied sites and empty sites. Such fully polarized states 
are appropriate only in the U — * oo limit, where double 
occupancy is fully suppressed by the Hubbard term and 
only the kinetic energy, -E k in = (Ht), remains relevant. 

In contrast to the spin case, where because of the SU(2) 
symmetry both the FM spin state and the AF spin state 
are unique, in the present orbital case without SU(2) 
symmetry there is already a plethora of physically dif- 
ferent ordered states even if one does not go beyond two 
sublattices. In particular, as shown by Takahashi and 
Shiba, 2 — Maezono and Nagaosa, 24 Shen et al.^ and par- 
ticularly stressed by Van den Brink and Khomskii^ it 
makes a big difference whether one builds an ordered 
state completely from complex orbitals (and empty sites) 
[i.e., -0, ipA, ipB = 0, 7r], leading to what we shall call com- 
plex states, or whether one uses exclusively real orbitals 
[i.e., ip, V'A; V\b = t/2], thus constructing real states. 
This can be conveniently demonstrated explicitly by for- 
malizing the description of the U — ► oo limit by means of 
the slave fermion formalism, which permits treatment of 
the general case (i.e., arbitrary ^i's and 0's). We present 
such states here in some detail, since the U — > oo limit 
will serve as a reference in the later discussion. 



So we introduce orbital bosons b 
to represent the occupied e 



t 
"I 

orbitals. 



(with rj 
\±)i = 



|0) = 



fey. I vac), and positively charged slave fermions fj to rep- 



resent the empty sites, |0), = //|vac). Thus the original 
electron operators are replaced according to c\ ± = b\ ± f i , 
and the Hamiltonian takes the form 



a (ij)\\a 



h_bl 



V "-■+"./ 
with the local constraint 



■^h_bl 



(3.3) 



(3.4) 



implementing the condition of no double occupancy. Or- 
bital order is then imposed by treating the bosons in 
mean field approximation, i.e. by making the replace- 
ments [compare Eq. (|2.6|) ]|2& 



b l+ -+ cos(^/2)e- ie '/ 2 , 
6i_ -> sin(^/2) e + 19 -/ 2 , 



(3.5) 



upon which the local pseudospin operators are given 
by Eq. (|2.7j> . We are then left with a Hamiltonian de- 
scribing fermionic holes moving in a background of fixed 
orbitals. 

In the case of FO order the result is explicitly 



njFO _ 
n U=oc — 



j*E E [l + 7sin^cos(0- X a)]/// 



a (ij)\\a 



(3.6) 

Upon Fourier transformation one obtains, reverting to 
the electron description, a single band with dispersion 
depending on the orbital angles {ip,9}, 



effl 0o (k) = -t A k +7sin^(cos0C k +sm0D k ) , (3.7) 



where 



G k 



cos k a + cos kb + cosk c , (3-8) 
i(2 cos k c — cos k a — cos kb), (3-9) 
i\/3(cos k a — cos kb)- (3.10) 



One notes that C k and Z? k transform as the 9 and e com- 
ponents of an E doublet, which makes Eq. (|3.7|) a cubic 
invariant (i.e., it does not change under the transforma- 
tion 9 — > 9 — 27r/3 and the simultaneous permutation 
k a — > kb — > k c — > k a ) . It will be useful to introduce also 



P k = ^2 e+lXa cos k a = C k + iD k , 



(3.11) 



which gets multiplied by the phasefactor e ~ l27T / 3 under 
the permutation k a — > kb — > k c — > k a , as well as the 
associated amplitude, 



# k = \Pk\ = {Cl + Dl}^ 2 = {cos 2 k a + cos 2 k b + cos 2 k c 

— (cosk a cos kb + cos kb cos k c + cosk c cos ka)} 1 ^ 2 , (3.12) 

which transforms as A\, i.e., has cubic symmetry^ 

Amongst the various phases with AO order let us consider first those of G-type (Neel- type), denoted by G-AO. One 
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obtains from Eq. (|2.3|) the two-sublattice Hamiltonian (i G A, j E B) 



-i/G-AO _ 
n U=oc — 



\t ^ ^ 1 1^ cos cos — i cos ip + sin 
+ 7[sin^+cos(0+-Xa) + « sin-0- sin(6> + -x Q ) \f}fj 



(3.13) 



depending on the orbital angles {ipA, ipB, 9a, Ob}, for which we introduce the shorthand notation for half the inter- 
sublattice angles, 



V>± = h(i>A±lpB), 



± 



(3.14) 



Upon Fourier transformation and diagonalization of the resulting 2x2 matrix this yields two electron bands (in the 
reduced Brillouin zone), 



-G-AO 
~-U=oo,± 



(k) = ±t| cos ip_ cos #-Ak + 7 sin-0 + ( cos #+Ck + sin 9 + D k 

cos ^ + sin 0_ Ak — 7 sin ?/>_ ^ sin 6* + Ck — cos 6* + L>k) j 



2-, 1/2 



(3.15) 



By a similar derivation one may obtain the electronic structure for the A-type and C-type AO phases. Using the 
same notation as above one finds 

£ c/=^°±( k ) = f j 1 + 7;j[sinip A cos{9 A - Xa) + sini/'s cos(6» B - x<xj] } cos/c Q 



a— a,fe 



^ cos ^_ cos + 7 sin ^ + cos + ^ cos ^4. sin — 7 sin sin cos 2 & c | 



(3.16) 



E 



± 



C-AO 
C/=oo,± 

7/ • 



(k) = t j ^1 + ^[sin^M cos6U + sink's cos^s]^ cosk c 



{ — (sin ?/>a cos 9a — sin ^ g cos 5 ) cos k c } ' 



+ | cosi/;_ cos 6L (cos /c a +cos fcfc) + 7sin?/; + ^cos (Q + ^) cos fc a + cos (9 + + cos k^j j 

/ 27T 27T \ 2l 1 / 2 ~| 

+ { cos V>+ sin 9- (cos fc a + cos fcj,) — 7sin-0- y sin (# + — J cos k a + sin (#+ + — ) cos fcbj } j 



(3.17) 

It is now straightforward to derive from Eqs. (|3.7p and ()3.15[) - (|3.17p the dispersion in any particular orbital-ordered 
phase with either complex or real orbitals. 

For the FOr [ip — tt/2] and the various AOr [ipA = 4>b = tt/2] rea/ states the dispersions are: 



4° r oc(k) = -*[A k +7(cos0C k +sin0D k )], 



_G-AO 
-C/=oo,± 



(k) 



_A-AO 
-C/=oo,± 



_C-AO 

-a=oo,± 



cos A k + 7 (cos # + Ck + sin # + D^j 
£(k) = -t{(2-7COs6L h cos6L)(A k -C k ) + 37sin0 + cos0_ £> k 
± ^7 sin 9 + sin 0_ (A k — Ck) + 37 cos 9 + sin 0_ £>k) 
+ ((cos6»_ +7Cos6» + )(A k + 2Ck)) 2 ] 1/2 }, 
(k) = ii{(l+7Cos(9 + cos6U)(/l k + 2C , k)± (7 sin 9+ sin 0_ {A k + 2C k 
2 cos 6_ - 7 cos (9+) (A k - C k ) + 37 sin 0+ Ac) | . 



(3.18) 
(3.19) 



(3.20) 



(3.21) 



In contrast to the complex states discussed below, all real states, whether FO or AO of any type and what- 
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ever the value of 9 (or 9 a and 9 b), explicitly break cu- 
bic symmetry, i.e., their dispersion is anisotropic. This 
nonequivalence between real and complex states is a man- 
ifestation of the broken SU(2) symmetry in the orbital 
Hubbard model (|2.9p . In extreme cases the dispersion is 
two-dimensional (2D). For instance, the dispersion of the 
"antiferro" (i.e. with Ta — —Tb) G-type AO state with 
alternating \x) and \z) orbitals (G-AOxz) [with 9a — 
and 6b = it] , 

efc»°t*( k ) = ±7«?k = ± 7 ^(cosfc a - cosk b ), 

(3.22) 

is 2D because the hopping along the c axis is fully sup- 
pressed when x 2 — y 2 and 3z 2 — r 2 orbitals alternate. 
Similarly, the dispersion of the fully |x)-polarized (FOx) 
state (9 — tt), 

e^ a 00 (k) = -t[A k - 7 C k ] 

= -t[(l+ ^j)(cosk a + cosfc 6 ) + (1 - 7) cos fc c ], (3.23) 

becomes 2D in the genuine orbital case (7 = 1), be- 
cause when only x 2 — y 2 orbitals are occupied, the only 
type of hopping allowed in this case, i.e. C-hopping [see 
Eq. (|2.10[) ], is suppressed along the c axis. 

Other states are also anisotropic, but typically have 
dispersion with contributions due to all three cubic di- 
rections. As an example, the dispersion of the G-type 
AO state with alternating 3x 2 — r 2 and 3y 2 — r 2 orbitals 
(6 A = — 9b = 27r/3) along the a and b cubic axes (G- 
AOab), 

e^f(k)=±t[-±A k + 7G k ] 

= =F^[(l+7)(cosA; a + cosfc b ) + (l - 27) cosfc c ], (3.24) 

is cubic at 7 = 0, but becomes predominantly but not 
fully 2D for 7=1. By contrast, the dispersion of the 



(9a = —9b = tt/2) state, with alternation between sym- 
metric and antisymmetric combinations, (|:r) + \z)) and 
(\x) - \z)), called G-AOsa, 

= ±7*[ — ^(cos k a + cos k b ) + cosfc c ] j (3.25) 

is quasi-one-dimensional (quasi-lD), qualitatively similar 
to that of the |.z)-polarized (FOz) state (9 = 0), 

4^(k) = -^ k+7 c k ] 

= -t[(l- 57) ( cos k a + cos k b ) + (1 +7) cosA: c ], (3.26) 

which becomes quasi-lD in the orbital case (7 = 1). 

The reduced symmetry of the FOx and FOz states is 
reflected in their respective densities of states, shown in 
Fig. [2fd), which lead to favorable kinetic energies [see 
Fig. EJd)], as discussed in Sec. IIII Bl Obviously such 
broken-symmetry states could be favored either in low 
dimensional systems, as the FOir state found for a 2D 
square lattice^ and suggested for bilayer manganites, 29 
or by a strong JT effect favoring a particular type of occu- 
pied e g orbitals due to oxygen distortions, as realized for 
instance in bilayer systems^ The latter applies also for 
the G-type AO states, which have typically smaller band- 
widths than the FO states; a few examples are shown in 
Figs. Mb) and He). 

As illustrative examples of the A-type and G-type 
phases with alternating real orbitals [either along the 
c axis or in the (a, b) planes], we give dispersions in 
each case for: (i) 9a = —9b = § , i.e. with alternat- 
ing (|a;) ± |z))/V2 states, and (ii) 9a = 0, 9b = n, with 
alternating \x) and \z) states, 



_A— AOsa 
~-U=oo,± 



(k) 



.A—AOxz 
'-U=oo,± 

_C-AOsa 
~-U=oo,± 

_C-AOxz 
'-U=oo,± 



(k) 
(k) 
(k) 



|t{2(A k - G k ) ± 7 [(A k + 2G k ) 2 + 

f r3 

cosfc a + cosfcb ± 7 

1 



3 

- (cos fc a — cos kb) 2 + cos 2 k c 



1/2 



'}■ 



^t(2 ± 7) (A k - C k ) = t(l ± i 7 ) (cos k a + cos k b ), 

k c ± ^7(cos fc a + cos k b ) 



i<[(l± 7 )A k + (2 T 7)G k 



= t 



COS fc c : 



A k + 2G k ± 7 [(A k + 2G k ) 2 + 9Dl] 1/2 } 

r 3 1 1 ' 

fc c ± 7 j-(cos fc a — cos kb) 2 + cos 2 fc c j 



(3.27) 
(3.28) 
(3.29) 

(3.30) 



The anisotropy of these phases is quite strong, and the Finally we consider the orbital ordered states with 
A-AOxz phase has even a 2D dispersion. complex orbitals. For two of these complex states, namely 
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the ferro |+)-polarized orbital order (FO+) [if) — 0] and 
the G-type alternating orbital order (AO±) with |+)/|— ) 
staggered orbitals [with ipA = and 4>b = tt], all cubic 
directions are equivalent, and one finds the dispersions 



_FCH 



,(k) = -tA k , 



and 



±jtB k , 



(3.31) 



(3.32) 



respectively. Thus, one finds that the dispersion of the 
FO+ state and its density of states, shown in Fig.[2ja), is 
that of a simple cubic lattice, as it originates entirely from 
the pseudospin-conserving hopping oc c\±Cj±, because 
at U — ► oo the alternating, non-pseudospin-conserving, 
hopping is fully suppressed by the imposed FO+ order. 
The reverse is true in the G-AO± state: here the disper- 
sion oc ±-Bk comes entirely from the alternating hopping 
the pseudospin-conserving hopping is fully 
suppressed by the AO± order. It is an important feature 
of both these complex states, built from cubic orbitals, 
that they retain cubic symmetry. Precisely for that rea- 
son these complex orbital ordered states were proposed 
as candidates for the ground state of the FM metallic 
phase of the manganites , 20 ' 23 to explain the observed cu- 
bic symmetry of the magnon spectra . 31 ' 32 

The other two orbital-ordered complex states break ex- 
plicitly cubic symmetry: the ^4-type and C-type AO± 
states. In the A-AO± state layers of |+) and |— ) orbitals 
alternate in the c direction, resulting in the dispersion 

4=i°i (k) = i* [(2 ± 7 )A k - 2(1 T 7)Oc] 

= t(cosfc a + cosfcb ± 7Cosfc c ). (3.33) 

This dispersion is qualitatively equivalent to that of the 
FO+ state, and thus the densities of states of the FO+ 
and A-AO± phases are the same. The reason is that 
replacing in the c-direction every second |+) orbital by a 
| — ) orbital does not affect the hopping parameter along c, 
and so the resulting doubling of the unit cell only halves 
the Brillouin zone without changing the dispersion. In 
the G-AO± state instead columns of |+) and |— ) orbitals 
alternate in the (a, b) planes, and one finds 



$=£±00 = it{A k + 2C k ± 7 



9Dl + {A k - C] 

1/2-, 



■)']} 

t [cos k c ± 7( cos 2 fc a + cos 2 kt — cos k a cos kb) 1 ^] (3.34) 

In contrast to the FO+, G-AO± and A-AO± phases, the 
G-AO± phase is not cubic symmetric. 

The densities of states of the complex states show 
a gradual crossover with increasing alternating orbital 
character from the full bandwidth of 6t for the FO+ and 
A-AO± phases, obtained also at U — both for the spin 
problem and for the e g band (Fig. [!}, to a narrower 
bandwidth of 2(1 + y/3)t for the G-AO± phase, and fi- 
nally to a bandwidth of At for the G-AO± phase. It is 
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FIG. 2: (Color online) Density of states N(u>) at 7 = 1 (in 
units of t — 1) for different orbital ordered phases: (a) com- 
plex orbital order: uniform FO+ [degenerate with A-AO±] 
(solid line), G-AO± (filled circles), and G-AO± (dashed line); 
(b) alternating real orbital order in G-type phases: G-AOsa 
(solid line) [the same density of states is obtained for C-AOsa 
phase], G-AOxz (dashed line), and G-AOab (circles); (c) al- 
ternating real orbital order in selected G- and ^4-type phases: 
C-AOxz (solid line), A-AOxz (dashed line), and A-AOsa 
(squares); (d) real uniform orbital order FOx (solid line) and 
FOz (dashed line). 



remarkable that, upon going from the FO+ phase to the 
A-AO± phase, the change from uniform to alternating or- 
bital order along only one cubic direction does not modify 
the density of states, while the density of states changes 
its shape completely upon going to the G-AO± phase, 
with a large spectral weight accumulated now close to 
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FIG. 3: (Color online) Kinetic energy gain (loss) AE/t with 
respect to FO+ phase [reference zero energy given by hor- 
izontal dotted lines], as a function of electron filling n for 
various orbital ordered phases (at U = 00): (a) complex or- 
bital order: G-AO± (filled circles), and C-AO± (dashed line); 
(b) alternating real orbital order in G-type phases: G-AOsa 
(solid line) [degenerate with C-AOsa phase], G-AOxz (dashed 
line), and G-AOab (circles); (c) alternating real orbital order 
in selected C- and A-type phases: C-AOxz (solid line), A- 
AOxz (dashed line), and A-AOsa (squares); (d) uniform real 
orbital order FOa; (solid line) and FOz (dashed line). The 
long-dashed line in (d) shows the kinetic energy for noninter- 
acting electrons in the e 9 band (disordered phase at U — 0). 



the band edges, resulting in a quite peculiar density of 
states with large maxima close to \u>\ ~ 2"/t, separated 
by a minimum with N(0) = at u = [Fig. g[a)]. The 
density of states for the C-AO± phase has a width of 



2(1 + \/3)t, and represents an intermediate case, having 
some features in common with that of the G-AO± phase. 



B. Densities of states and kinetic energies in 
orbital ordered states 



It is worthwile to consider next the densities of states 
of various orbital ordered states in a little more detail 
(see Fig. [2]), and investigate their consequences for the 
kinetic energy (Fig. [3]) . Focusing first on the bandwidth, 
we note that for any FO state this takes the maximum 
attainable value 6t. This result is not limited to the FO 
states considered explicitly above, i.e., the complex FO+ 
and the real FOz and FOa: states, for which it was already 
pointed out by Van den Brink and Khomskii, 20 but holds 
in general, i.e., for arbitrary ip and 9, as readily shown 
from Eq. (|3.7[) . Moreover, the result is independent of 7 
and so holds both in the orbital case and in the spin case. 

By contrast, in any G-type AO state the bandwidth is 
smaller than Qt [see Fig. dfb)], and depends on 7. In 
particular, in any G-type "antiferro" state (with Ta = 
—Tb, so ipB — 7T — ipA and 0b — ®a — tt), such as the 
complex G-AO± state or the real G-AOxz and G-AOsa 
states considered above, the width is proportional to 7 
(viz. 4jt, 2v / 37t, and A'yt, respectively for those three) 
as follows from Eq. (|3.15p . In such a state the bandwidth 
therefore, correctly, collapses to zero in the spin case (7 = 
0) where hopping is completely suppressed by the AF 
spin order^ The important point to note here is that 
at finite 7, and so in particular in the genuine orbital 
case, the bandwidth even of an "antiferro" state is finite 
though smaller than that of the FO states. Thus, while 
in the spin case the kinetic energy of carriers is fully lost 
when going from FM to Neel-type AF order, this is not 
so for the analogous FO to G-type AO transition in the 
orbital case. 

One might still be tempted to believe that, as fa- 
miliar from the spin case, also in the orbital case FO 
order is most favorable for lowering the kinetic energy 
of charge carriers, simply because this gives the largest 
bandwidth. However, the situation is not that simple, 
not only because there are several inequivalent FO states 
with different densities of states which have nevertheless 
the same bandwidth, but also because some A-AO and 
C-AO phases have again the same bandwidth, and so one 
really has to consider the details of the density of states 
in each case. This is demonstrated in Fig. [3l which shows 
the kinetic energy gain AE with respect to the complex 
FO+ state as a function of electron filling n for various 
FO and AO states with (complex or real) orbital order, 
obtained by straightforward integration of the respective 
density of states. Indeed, at small electron filling n, and 
also at small doping x = 1 — n, AE is lower for the (FO, 
A-AO and C-AO) states with full bandwidth 6t than for 
any state with a narrower density of states, in particu- 
lar for the G-type AO states of Fig. E^b), because the 
first doped holes enter in the former case with an energy 
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~ — 3t close to the band edge, while the lowest accessi- 
ble energy is higher in all G-AO states. Note that the 
orbital order observed in LaMnC>3 is close to that of the 
C-AOsa phased and this phase has the same density of 
states as the G-AOsa phase [see Fig. Efb)], and thus has 
a rather unfavourable kinetic energy [Fig. E^b)]. This 
demonstrates that both an interplay between spin and 
orbital order due to the SE interactions at finite U, and 
the JT interactions between orbitals on neighboring sites, 
induced by the coupling to the lattice, play an important 
role in real materials and stabilize the orbital order ob- 
served in undoped LaMnOs&S^. 

Among the states with AO order of real orbitals, but 
FO order along one or two cubic directions, we identi- 
fied three phases, C-AOxz, A-AOxz, and A-AOsa, which 
have lower energies than the FO+ phase close to n = 
and n = 1 [Fig. 02c)]. All of them have the full bandwidth 
6t [Fig. [2Jc)], but a finite density of states at uj = — 3t 
gives the A-AOxz phase the lowest energy of these phases 
at very low n or x. At somewhat higher filling n ~ 0.07 
(doping x ~ 0.07) the other two phases take over, and 
are in fact more stable than the FO+ phase in the en- 
tire regime of n. This follows from the large densities of 
states of these phases at \oj\ ~ It. In contrast, the A- 
AOxz phase with a large spectral weight close to uj = 
has a higher energy than the FO+ phase in a broad range 
of electron filling 0.11 < n < 0.89. 

The above discussion shows that at finite but still 
rather modest electron filling or doping, the overall shape 
of the density of states becomes more important, and the 
states with large density of states near the band edges 
could be favored a priori, even in cases when the band- 
width is smaller that 6t. An interesting example here 
is the complex "antiferro" G-AO± state, with its en- 
ergy falling below that of the complex FO+ state for 
n > 0.27 or x > 0.27 [Fig. [5fa)], because of the large 
number of states available in the G-AO± state just close 
to the band edges at \u>\ = 2t, whereas in the FO+ state 
the energy of available electron states, though initially 
— 3t, rises rapidly with increasing doping [Fig. 02a)]. 
However, in reality the transition from FO+ to G-AO± 
state does not happen, as the real FO (FOz and FOx) 
states have even lower kinetic energy throughout than 
both complex states. This can be ascribed to the lower- 
dimensional nature of their dispersion and the resulting 
different location of the Van Hove singularities, which 
[compare Fig. 03d)] enhances the density of states near 
the band edges at ±3i and at the band center for the 2D 
FOx state, and in the intermediate range t < \u>\ < 2t 
for the quasi-lD FOz state. As a result, at small filling 
(doping) the kinetic energy gain AE is the lowest one for 
the FOx state, while at larger filling n > 0.30 (doping 
x > 0.30), the FOz state takes over. However, in this 
regime of electron filling the energy gain AE for the A- 



AOsa phase is lower by a few percent, and the two phases 
may be considered as practically degenerate. 

For comparison and later reference we have included 
in Fig. 03 also the kinetic energy for the uncorrelated e g 
band (the correlated OL phase is analyzed in Sec. 021- Of 
course, at U = any kind of orbital order is absent and 
one finds by far the lowest kinetic energy for the disor- 
dered e g orbitals. The e g bands have then the dispersion 
given by 

£c/=o,±(k) = -t(A k ± 7.B k ). (3.35) 

Remarkably, these bands at U — represent formally a 
superposition of the FO+ and G-AO± bands at U = 00, 

s U= o,± (k) = 4°+ (k) + E%zi°i (k). (3.36) 

and so naturally also show full cubic symmetry and a 
bandwidth equal to 6t [Fig. Olid)]. One notes that, be- 
cause both pseudospin-conserving and non-pscudospin- 
conserving hopping channels fully contribute here, con- 
siderably more kinetic energy can be gained than in any 
of the orbital-ordered states. In particular, as Fig. OUd) 
shows, there is a large density of states at and near the 
band edges, and thus AE is the lowest in this disordered 
state already at small electron filling n, and then remains 
so throughout. Of course, this large kinetic energy gain 
will be partly lost for large U near n = 1, where at least 
one hopping channel gets partially suppressed by electron 
correlations. However, the result here indicates that the 
tendency towards the OL state with disordered e g or- 
bitals is particularly pronounced. We shall come back 
to this point, presenting more evidence in favor of the 
correlated OL phase, in Sec. 021 



IV. HARTREE-FOCK APPROXIMATION 

A. Instabilities towards orbital order 

We turn now to the orbital Hubbard model (|2.9p with 
finite U, where it is to be expected that polarization, 
when it occurs, need not be complete but can be par- 
tial, as in the spin case. Also, the existence of orbital 
ordered states will in general require a sufficiently large 
U/t. Which instabilities towards orbital ordering oc- 
cur and at what value of U jt can be investigated either 
by considering the corresponding susceptibilities, e.g. in 
random phase approximation^ or by comparing the en- 
ergies determined in the HF approximation^ If various 
ordered states are possible, one needs to calculate their 
energy (or free energy at finite temperature) to determine 
which one is actually realized. 



In the absence of SU(2) symmetry it is not sufficient to decouple the interaction term in Eq. (|2.9|) in the familiar 
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mean-field way, nj+n,- ~ ((n,-_(-)ni_ + n,+ (nj_) — (ni+)(nj_)), but one needs instead the general HF decoupling, 

((n l+ K_ + m+im-) - (n i+ )(n < _)) - ((^ + )c|_c i+ + 4 + q_(^r> - W<IT>) • (4-1) 



In the FO case, i.e., when one assumes a single three-component order parameter, T z — (Tf ), T + = (T+), TL = (Tj ), 
one obtains upon Fourier transformation the HF Hamiltonian 



FO 
HF 



= E (4 + 4 



Wn-UT z -tA v 



-UT_ - -ftP* 



k v -UT+- 1 tP^ iUn + UT z -tA u 

~\Un 2 + U(T 2 + T+T_), 




with P k given by Eq. (|3~PTj) . The eigenvalues are (with T+ = Te ie , T_ = Te~ l9 , so that T+T_ = T 2 = T 2 + T 2 ) 



4°(k) = -^ k + C/(in±£; k ), 



where 



T 2 + T 2 



2^T(cos0C k + sin<9D k ) + (^)~ B ] 2 



1/2 



and the HF groundstate energy per site is then given by 



E ™ = ^E [Mk)e ro (k) +n + (k)4°(k) 



--[/n 2 + [/(T 2 +T 2 ), 



r 



(4.2) 

) 

(4.3) 
(4.4) 

(4.5) 



where n_(k) [n + (k)] is the occupation number of the 
lower (upper) band. For large U (> 6t) a gap opens, 
and so for less than half-filling only the lower band is 
occupied. Setting the derivatives of E^p with respect 
to n, T z , T, and 9 equal to zero then yields the self- 
consistency equations 



^E-(k), 



T z = 



N 

k 

2N^ n 1 V 



T 



= ily n _( k ) 



2N 



— V n_(k) 



T+ 2l( cos 6»G k + sin (9 L> k ) 

T(sin0C k -cos0£> k ) 
Sk ' 



(4.6) 
(4.7) 

(4.8) 
(4.9) 



While Eq. (|4.6p is trivially satisfied in the sense that it 
simply fixes the Fermi level for given filling n, two general 
conclusions can be proven from the remaining three equa- 
tions. Firstly, it follows from Eq. (|4.9p , because of the 



dependence of i? k on C k and D k [see Eq. (|4.4p ] and the 
explicit form of the latter two functions [see Eqs. (|3.9|) 
and (|3.10p ], that for nonzero T the azimuth must equal 
either (or equivalently +27r/3 or — 27r/3) or n (or equiv- 
alently — ir/3 or +7r/3), i.e. the projection of the order 
parameter on the 'real' equatorial plane has to be along 
one of the cubic directions. Secondly, it follows that both 
a purely complex state (i.e., T — 0, T z ^ 0) and a purely 
real state (i.e., T z — 0, T ^ 0) are permissible states, 
in the sense that T z = is a self-consistent solution of 
Eq. P~7|) and alternatively T = is one of Eq. ([3~8"|) . We 
remark that both these properties of the possible states 
need not be postulated or assumed but are proven here 
from the HF self-consistency equations. 

As C-type and ^4-type AO phases would give qualita- 
tively similar results, we will consider from now on only 
G-type AO phases, and denote them for brevity by "AO" 
instead of by "G-AO" . So we assume independent three- 
component order parameters on interlacing sublattices A 
and B, T Z A = (Tf) A , Tf = (Tf) B , etc. Then the HF 
Hamiltonian is 



(J 



-Hi 



AO 
HF 



E 

k 



^A,k+ 
,t 

-A.k- 

L B,k+ 

J 



( \ Un 



\ 



\Un 2 



Ll.k- / 



-UT£ 
~tA k 



UT A 

7^k 
-tA k 



-iA k 
\Un- UTf 



-tA k 

_ UT B 



(T?) 2 +T*T*]. 



"A,k-| 



c A,k- 
C B,k+ 
V C B , k - J 



(4.10) 
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Like above, the HF groundstate energy per site is then formally given (with T A = T A e l6A , etc.) by 

E ™ = W EE^W^W " 3^ 2 + M(T A ) 2 + (T 2 A ) 2 + (T B ) 2 + (TB) 2 ] , (4.11) 

13=1 k 



where the sum on (3 is over the four bands and that on k 
is over the reduced Brillouin zone. However, as the 4x4 
matrix in Eq. (|4.1ip cannot be diagonalized analytically 
in the general case (i.e. for arbitrary order parameters), 
no further progress can be made like in the FO case. In 
particular one cannot strictly prove that purely real or 
purely complex states are permissible solutions. 

Yet this still seems likely, and if one makes this as- 
sumption, then for the case of the complex (AOc) state, 
i.e. with T A = T B = 0, the 4x4 matrix simplifies enough 
to obtain explicit expressions for the band dispersions, 



r AOc 

£ /9 



(k) 



-U 



where 



F 2 ,k = 



T B \ 2 



I 2 

- d k 



1/2 

1 

1/2 



(4.12) 

(4.13) 
(4.14) 



Setting the derivatives of E^p ° with respect to n, T A and 
T B to zero yields again HF self-consistency equations. 
From these one easily proves that T A — —T B , i.e. that 
the stable complex state is actually the AO± state. 

For the case of a real (AOr) state, i.e., with T A = 
Tf = 0, an analytic solution is also possible, but this is 
so unwieldy as to be impractical. However, if one further 
assumes that T A = T B = T one can derive the approxi- 
mate expressions 



e AOr (k) = ±t[j4 k cos0_ + 7 (C k cos0 4 



U(±n±G k ) 



Ac sin 0+)] 
(4.15) 



where 



T 2 + 









'{Al 



sm 0- 



-7 2 (C k sin6» + -D k cos6»+) 2 } 



1/2 



(4.16) 



valid in the large U limit (U/t 3> 1), and again obtain 
analytic self-consistency equations by taking the deriva- 
tives of E AOr with respect to n, T, 0+, and 6L. From the 
latter two one can now prove the following. First, that 
9 + = 0, i.e., 9a = —&B-, so the pseudospin vectors on the 
two sublattices are mirror images of one another with 
respect to the cubic direction 8 — (or the equivalent 
ones 9 = ±27r/3). Second, that cos6L ~ —(U/3t)x for 
x <C t/U, i.e., 9a — 7r/2 + (U /3t)x, so that at zero doping 
the stable solution is the AOsa state, and with increasing 
doping the pseudospin vectors tilt slightly away from the 
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FIG. 4: (Color online) Stoner instability towards FO+ partly 
polarized states (full line) as a function of band filling n in 
the orbital Hubbard model (7 = 1), and the inverse of the e g 
density of states of Fig. 0Jd) (dotted line); full polarization 
occurs only in the limit U = 00. For the spin model (7 = 0), 
the corresponding Stoner instability (dashed-dotted line) is 
given by the inverse of the density of states shown in Fig. 
(Ha), while saturated FM states occur above the dashed line. 



cubic direction, making the solution gradually resemble 
more the AOab stated 

As an example of the HF instability at intermediate 
U we have investigated how the complex FO+ (or the 
equivalent FO— ) state develops when U increases, using 
Eqs. (|4.2|) and (|4.5[) . First, at 7 = one recovers the 
Stoner criterion U N{Ep) — 1 for the onset of the FM 
order with increasing U, with the FM saturated states 
becoming stable at still larger but finite value of U (Fig. 
[4]). By contrast, in the orbital model at 7 = 1 the insta- 
bility is qualitatively different, and the FO+ (FO— ) state 
appears as a global property of the band rather than as an 
instability at the Fermi surface. The instability occurs at 
higher values of U for any filling than in the spin case — 
actually the value of the critical U is very close to that 
giving full magnetic polarization in the spin case. 

Here, unlike in the spin case, the FO order implies 
that the electronic bands are changed — they develop 
an additional splitting above a critical value of U, which 
modifies the shape of the bands and leads to a finite 
order parameter T z = (Tf) ^ 0. This mechanism of 
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the instability resembles that known in the spin case for 
the onset of antiferromagnetism. The critical value of U 
above which weak order appears has therefore no relation 
to the actual shape of the density of states (see Fig. 

We decided not to investigate the phase diagram of 
the orbital Hubbard model in the HF approximation in 
detail. Instead, we concentrate first on the qualitatively 
novel aspects of various possible ordered states in the 
regime of large U , where, as we will see, the contrast 
with the spin case manifests itself in the most trans- 
parent way. Using these results, we will then comment 
of the HF phase diagrams analyzed in detail by several 
groups, 



^24^ i n Sec. EVD] 



B. Superexchange in the complex orbital states 

We have already seen that the analysis of the orbital- 
ordered states simplifies when the splitting of the quasi- 
particle bands cx UT Z or cx UT is sufficiently large that 
it opens up a gap and only the lower band (lowest two 
bands for AO order) is (are) partly occupied when n < 1 . 
It is then straightforward to calculate the energy and the 
order parameter by summing over the occupied states. 

Consider first the ordered states with complex orbitals. 
In the case of the FO+ state the equation for the order 
parameter, from Eqs. (|4.8[) and (|4.4|) . takes the simple 
form (because only T z 0, while T = 0): 



O N ^ 



n_(k) 



2N ^ E k 

k 



£k 



1 



7* A 2 R 2l 1/2 

-Dl. 



UT Z 



(4.17) 
(4.18) 



Equation (|4. 18|) shows explicitly that, unlike in the spin 
case, T z = n/2 only at U = oo, basically because the 
saturated FO+ state is not an eigenstate of the orbital 
Hubbard model given by Eq. (|2~9)) . Thus the FO+ state 
is again seen to resemble the AF phase in the spin model. 

Similarly, in the AO± phase for large enough U the 
order parameter is given by 



= 1 l_^pn 1 (k)+n 2 (k) 



1 



UT Z 



At 



1/2 



(4.19) 
(4.20) 



rather similar to the FO+ case (|4.18|) . but with the inter- 
change 4 k <-> 7i?i c . The reason is readily recognized from 
Eq. ()2.9|1 : for FO+ order, the diagonal hopping cx cJ ± c J± 
that gives A^, is order- preserving, while the off-diagonal 
terms cx c| ± c J=F that produce £?k are order-perturbing 
and reduce T z . For AO± order this is reversed: the off- 
diagonal hopping cx 7 that gives £?k is compatible with 
the order, while the diagonal one that gives A^ disturbs 
it. 

The similarity between the FO+ and AO± states at 
7 ~ 1 becomes even more transparent at large U (i.e. 



^ t), where near half- filling (i.e. for small x — 1 — n > 0), 
upon expansion up to first order in t/U, 



rpFO+ 
rpAO± 



3 7 2 



1-x) 



(1-x) 2 
3 - 2x 
(l-x) 2 



Note that a SE contribution cx ("ft) 2 j\J appears also in 
the FO+ state, because the off-diagonal hopping permits 
virtual charge fluctuations. This result is again quali- 
tatively different from the spin case, where the SE con- 
tributes only in the AF states, and so destabilizes uniform 
FM spin order. In the genuine orbital case (7 = 1) the 
reduction of the order parameter by SE is the same for 
FO+ and AO± at x = 0, but at x > it is slightly larger 
for the FO+ phase. The corresponding expressions for 
the energy, up to second order in t/U , become 



^ FO+ = -^E"-( k )^ 



37 2 



t 



E 



AO± 



2(1 -x) U 

-^EM k )- n 2(k)]i?k 

k 

3 - 2x t 2 
2(1 -x) U' 



77- ( 4 - 23 ) 



(4.24) 



Both are seen to be composed of the U — 00 kinetic en- 
ergy (compare Eqs. (|3.31[) and (|3.32p for the dispersions) 
and a (negative) SE energy. Surprisingly, near half-filling 
the energy per site of the FO phase is lower than that 
of the AO phase at any value of U, not only because the 
FO phase gains more kinetic energy cx — 3tx than the AO 
phase cx — 2tx, but also because it has lower SE energy. 
Instead, AO± order yields lower energy at larger doping 
x > 0.27 as a consequence of its peculiar density of states 
[Fig. ^a)] i? '- i — Note that this is opposite to the spin case 
(7 = 0), where the Neel (AF) state has lower energy near 
ri = 1 and the FM state takes over only above a critical 
doping x c ~ t/2U. 

We emphasize that we have compared as yet only the 
two complex states with one another, with the express 
purpose of contrasting the behavior of these orbital states 
with that of the corresponding spin states. To establish 
what the most stable orbital-ordered state is, we still have 
to consider the real states. 



C. Superexchange in the real orbital states 

The results obtained for the ordered phases with real 
orbitals are qualitatively similar. We focus here on the 
representative cases of the FOx, the FOz, and the (G- 
type) AOsa states, which we have shown in Section IIVI 
to be solutions of the HF equations. Note that the AOsa 
phase is representative for G-type AO order. For sim- 
plicity we ignore here the small higher order correction 
to the equations below j^l which occur when the actual 
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occupied orbitals deviate from those of the AOsa state 
towards those pertaining to the AOab state as discussed 
above. 

At large U/t one finds near half- filling for the order 
parameters 

i 




e fo z = - t ±y n 

k 
k 



The corresponding energies in these ordered phases are 



2 j-2 



— X 



_(k)(^ k - 7 C k ) 

_Ck)(A k + 7 c k ) 

Wk)-n 2 (k)1 7 C k - 



4(1 - x) U' 

3 7 2 t 2 
4(1 -x) U' 

6-4x + 3 7 2 t 2 



(4.27) 
(4.28) 
(4.29) 



Unlike the complex states, the real states are seen not 
to be degenerate in the undoped case x = 0. The AOsa 
state has the lowest energy here, even though the SE 
contributes also in the FO states. However, we find the 
same qualitative difference with the familiar AF and FM 
states for spin order as we found for the complex orbital 
states — again the SE contributes both in FO and in AO 
states. 



Finally, we remark that the SE contributes also in any 
other phase, either with mixed FO and AO order (e.g. 
in the C-AO and A- AO phases of Sec. IHIj) . or in a dis- 
ordered OL state. Depending on whether the occupied 
orbitals on a given bond are identical or not, virtual pro- 
cesses due to pseudospin non-conserving or pseudospin 
conserving hopping contribute, and we have verified that 
qualitatively similar results are then obtained to those 
presented in Eqs. (|4.25M.29|) above. Such terms would 
play a role in the low-doping regime and would deserve 
a separate study in order to establish the phase diagram 
of weakly doped manganites. Note that in that regime 
also the spin-dependent SE plays a prominent role, and 
the present orbital Hubbard model ()2.9|1 , which implicitly 
assumes FM order, becomes insufficient to describe the 
physical properties of the real materials. On the other 
hand, the SE terms, being all oc t 2 /U , vanish in the limit 
of large U which we consider in Sec. V, and hence they 
have no consequences for the stability of the OL phase 
at U = oo. 



D. Qualitative understanding of the Hartree-Fock 
phase diagram 

Finally, let us analyze the possible instabilities of the 
orbital Hubbard model (|2.9p in the HF approximation. In 
the large U limit relevant for such instabilities, the total 
energy consists of the kinetic energy at U — oo, discussed 
in Sec. IIII Al and a negative SE energy. While we do not 
intend to make a quantitative comparison between the 
various phases stable in the HF approximation, knowing 
that they are anyway destabilized by the correlation ef- 
fects (see Sec. |V}, this now enables us to get a simple 
interpretation of the HF phase diagram of the genuine 
e g orbital model (7 = 1).2&2 & 24 i 25 using the large U ex- 
pansion. These earlier HF studies have shown that at 
half-filling, and in the regime of small doping, for U > 6t 
the most stable state is the real 'antiferro' orbital state, 
with the orbitals close to those found in the AOsa phase. 
In this regime the SE energy dominates, and indeed the 
largest energy gain is then given by Eq. (|4.29[) . At in- 
creasing hole doping, however, the kinetic energy of holes 
moving in the FOx background is much lower than that 
in the AOsa phase (see Fig. [3]), leading to a transition 
to 'ferro' orbital states when the difference between the 
SE terms (xt 2 /U is overcome by the difference between 
the kinetic energies of these two phases. The region of 
the AOsa phase in the phase diagram decreases when 
the SE gradually looses its importance with increasing 
U, as shown by the numerical result of Van den Brink 
and Khomskii. 20 

At U = 00 the FOr order is found in the HF approxi- 
mation at any doping x > 0. However, at large but finite 
U the SE is larger in the FO+ than in either FOa; or 
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FOz phase, while the difference in the kinetic energy is 
small [Fig. Eld)], and thus the FO+ state is the first sta- 
ble 'ferro' state at intermediate values of 8 < U/t < 12 
and x ~ 0.15. However, when x increases further, the ki- 
netic energy difference between the FOx and FO+ phase 
dominates, and the orbital order changes to FOx. As 
the SE energy of the two real FOx and YOz states [see 
Eqs. (|4.27p and (|4.28[1 ] is the same, the difference in the 
kinetic energy gives a second transition from the FOx to 
the YOz phase with increasing x. At small and interme- 
diate U/t < 12 one finds eventually at x ~ 0.5 the AO± 
phased which is stabilized in this regime by a combined 
effect of large SE energy gain and low kinetic energy (see 
Fig. [3|) which follows from the peculiar density of states 
of this phase. 

In a 2D model the phase diagram is quite different,— 
and is dominated by the generic tendency towards x 2 — y 2 
polarization within an (a, b) planed The AO order is 
then followed by the FOx phase above a critical dop- 
ing, which decreases with increasing U/t. We note that 
the region of the FOx phase is enlarged by the offdiago- 
nal hopping terms cx 7^,— in agreement with the above 
observation that these terms stabilize the FO phases at 
finite U due to the respective SE energy contributions. 



V. ORBITAL LIQUID STATE 
A. Kotliar-Ruckenstein slave boson representation 

To understand further the essential differences between 
orbital and spin physics, we develop now an approximate 
description of the correlated OL disordered state. This is 
of crucial importance as the HF approximation permits 
only a comparison of ordered states with one another, 
and therefore does not allow to draw any conclusions con- 
cerning the stability of the orbital-ordered states with re- 
spect to disordered states. This is well known from spin 
models — for instance, the FM states in the 2D Hub- 
bard model are stable only in a narrow range of doping 
x < 0.29 near half-filling, 39 while the HF approximation 
predicts FM to be stable at any electron filling n. 

We will argue below that indeed orbital (FO or AO) 
order is not robust at 7 = 1 and gets replaced by a 
disordered (OL) phase, if one goes beyond the HF ap- 
proximation and includes electron correlation effects. As 
we have already seen, the orbital problem is richer than 
the spin case, as various ordered states are nonequiva- 
lent when the SU(2) symmetry is absent. Therefore, we 
shall consider only the limit of very strong correlations 
and investigate the stability of orbital order specifically 
in the U = 00 limit, where the OL competes with fully 
saturated FO [see Eqs. P~2l"j) and (|3~2"5) ] and AO [see 
Eqs. P~2"2"|) and P~2"5] states. 

In order to obtain a reliable variational method to cal- 
culate the correlation energy, we have followed the slave 
boson approach introduced by Kotliar and Ruckenstein 14 
for the spin Hubbard model, and have adapted it to the 



orbital case. In this approach the Fock space is enlarged 
by the introduction of three auxiliary bosons at each site, 
one for each local configuration, viz. bi + and asso- 
ciated with the single-occupancy configurations \i+) and 
\i— ), and e% with the empty configuration \i0) (double 
occupancy is excluded at U = 00). Then a physical 
fermion (electron) c is represented by a pseudofermion 
/ and two accompanying bosons according to an expres- 
sion like ctg = //a;&Jge,-, where the two bosons keep track 
of the change of the local configuration when an electron 
is added. 1 This construction, however, must preserve the 
cubic symmetry of the Hamiltonian (|2.9p . implying that 
it has to be gauge invariant with respect to those U(l) 
rotations in orbital space that correspond to a permuta- 
tion of the cubic axes. The relevant rotation operator is, 
for arbitrary rotation angle 9, 

11,(6) = exp(-i0Tf). (5.1) 

The complex orbitals pick up just a phase factor under 
any rotation of this form, and the operators {cL_, cj_} 
transform as 

umcl+uho) = e-^ 2 4 + , 

Um4-U}{e) = e+ i$ / 2 4_. (5.2) 

As already indicated in Section |TT1 the orbital Hubbard 
Hamiltonian (|2.9[) is invariant under a uniform rotation 
at all sites, if the common rotation angle 9 is one of the 
three cubic angles —4tt/3, +4tt/3, 0, and if this is accom- 
panied by a corresponding shift of the "gauge angles" Xa 
by +27r/3, — 27r/3, 0, respectively. Actually the diagonal 
hopping terms in (|2.9[) are invariant under the U(l) trans- 
formation (|5.ip even for arbitrary 9, as a consequence of 
the SU(2) symmetry of the spin Hubbard model, while 
the off-diagonal hopping terms pick up phase factors, 

c i+ Cj_ h-> e l9 c\ + c j _, 

4-S-+ " e +ie 4_c J+ , (5.3) 

which get compensated by the shift of the Xa if is a 
cubic angle. As the three cubic-angle transformations 
amount to a forward and to a backward simultaneous 
cyclic permutation of axes and orbitals and to the iden- 
tity, respectively, the invariance expresses the cubic sym- 
metry of the Hamiltonian. 

Therefore, we take the slave boson representation as 

4 = b\ ± f^ ei , (5.4) 
corresponding to a representation of the local states by 

1*0) = et|vac), 
\i+) = 4 + \i0) = &t + /^|vac>, 
\i-) = 4_\i0) = b\_fl\vac), (5.5) 
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and we impose that the boson and pseudofermion oper- 
ators transform under U(l) rotations^ as 



Ui(0)bt ± u}(9) = e^b. 



i± i 

umfLuHe) 



(5.6) 



Note that the phase of the boson operators b\ ± changes 
twice as fast as the phase of the pseudofermion operators 
fj ± , i.e. the bosons have pseudospin T = 1, while the 
(pseudo)fermions belong to T — 1/2. This property guar- 
antees that the U(l) rotation behavior of the electron op- 
erators, as given in Eqs. (|5.2|) . is correctly reproduced by 
the transformation (|5.4p . Thus the present formulation is 
indeed gauge invariant and preserves the cubic symmetry 
of the orbital problem, like the SU(2)-invariant formula- 
tion introduced by Fresard and Wolfle preserves the full 
rotational symmetry for the spin system^ Clearly, the 
construction of a gauge invariant formulation is greatly 
facilitated by our use of the complex-orbital represen- 
tation, but a similarly gauge invariant representation in 
terms of real operators can also be constructed, and is 
given in the Appendix. 

The enlarged Fock space contains also unphysical 
states which must be eliminated by imposing con- 
straints as in the original formulation by Kotliar and 
RuckensteinJ^ 



bib 



i+ 



-b\_b t _+4e t = 1, 



b\ + b i+ = f\_f i _, 6l fo, = flfa 



(5.7) 



and implemented by means of Lagrange multiplyers 
{Xi, Mi-}- The first constraint excludes double oc- 
cupancy, the other two eliminate the unphysical singly- 
occupied states b\,j |vac) and | vac). The elec- 
tron density and the z-component of the pseudospin can 
then be described at each site either by slave boson or by 
pseudofermion operators, 



uncorrelated (U — 0) limit. The renormalized boson fac- 
tors take the form 



sJ(l-4^-b\ T b lT )(l-bi ± b i± ) 



(5.12) 



where it is important that the operator expression under 
the square root in the denominator is U(l) invariant, so 
that z\, (zj_) transforms under (|5.6|) exactly as b i+ (b\_). 
Then the Hamiltonian in the slave boson representation 
at U — oo becomes 

n u=™ = E [4+fl-fj- z j+ + 4-fl+fj+ z i- 

a (ij}\\a 

- ^ (&J + & i+ + + eje, - 1 



iX 



with A = ± and A = — A. The Hamiltonian commutes 
with the constraints and thus does not connect the phys- 
ical and the unphysical subspaces of Fock space. 

In the mean-field approximation we replace the boson 
operators by their averages. In order not to spoil the 
cubic invariance only their amplitudes are replaced by 
c-numbers, while their phases are prescribed to behave 
still according to Eq. (|5.6|) . 43 So we set for the boson 
invariants 



(bl+k, 



bl, 



(b'M 



(4 e i) = e ?. 



(5.14) 



where and e$ are real quantities, i.e. do not 

contain any nontrivial phased For the offdiagonal, non- 
invariant, two-boson products we set 



b\ + h + 



fi+ fi+ + fi—fi—i 



(5.8) 



Wi + bi + -b\-b i _) = \{flf i _-fl + h + ). (5.9) 



The other two components of the pseudospin operator 
can only be represented as 



T; 



bl+k-fLf^ 



T- =b]_b i+ flf i _, 



(5.10) 
(5.11) 



and cannot be reduced to expressions in terms of either 
slave bosons or pseudofermions alone.— 

As in the spin case one further has to renormalize the 
bosonic factor in Eq. (|5.4|) in order to recover, when a 
mean-field approximation is going to be made and the 
constraints are no longer rigorously obeyed, the correct 
unrenormalized hopping for the pseudofermions in the 



e~ iSi , (5.15) 



{b\_e i )=b i .e i e+ i - 9 \ (efo.) = h 



where the 'phase operator' is understood to transform 
as 



Ui(0) diU}(6) = di + 6, 



(5.16) 



and in particular assumes the cubic values t? a , $6, and 
f9 c when the two-boson operator product occurs in an 
expression taken along the a-axis, 6-axis, or c-axis, re- 
spectively. The last average of Eqs. (15. 14|) controls the 
number of holes in the e g band, ef = x, for a phase with 
uniform charge density. The constraints give then the 
following self-consistency conditions, 



bl 



(flfi-h 



bl = (fj + f i+ ), 

l-x, (5.17) 
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while the renormalization factors become 



4) 



with 



1 - (fW 



1 - (n 



(5.18) 
(5.19) 



The exponentials containing can be eliminated from 
the Hamiltonian by absorbing them in the pseudo- 
fcrmions, according to 



ft 



(5.20) 



Note that this definition ensures that the f}± transform 
properly under U(l) in accordance with Eq. (|5.6[) . 

Within the slave boson mean-field approximation one 
thus finds an effective Hamiltonian for pseudofermions 
subject to local constraints, and with renormalized hop- 
ping. In the case of orbital-ordered phases its precise 
form depends on the assumed type of state, with the 
hopping renormalization factors q i± either uniform or al- 
ternating between two sublattices. Here we present only 
its simpler form, adequate for uniform phases, such as 
FO and OL states, in which the renormalization factors 
and Lagrange parameters can be taken site independent, 

1 M«+ 

<y>l|a 



n 



MF 



\ " 



2 ^ 



q-fll- 



j+ 



(5.21) 



with fii\ = f\\fi\- The present formalism reproduces 
the results of Kotliar and Ruckenstein for the spin model 
(7 = 0) with hopping ^i, and gives the same results as 
the Gutzwiller approximation^ and so q + and g_ will 
be called also Gutzwiller factors. 

The ordered states can be obtained within the present 
KR slave boson approach by a proper choice of the La- 
grange multipliers. For instance, the FO+ state is now 
obtained from Eq. (|5.21| by imposing (fii-) = by 
means of the condition = —00 (while /i + = 0). Such 
states do not experience any band narrowing, as dou- 
ble occupancy is rigorously eliminated at U - co, and 
the correlation energy vanishes^ As a result, only the 



-FO 



- — tAk band is partly filled in the FO+ state, 
while the eff^^ ± (k) = ±jtB^ bands are filled in the 
AO± state. Real orbital-ordered states can also be ob- 
tained, using the formalism described in the Appendix. 
Therefore, in the U = 00 limit one reproduces the re- 
sults of the HF approximation described for these states 
in Sec. llVl 



B. Nature of the orbital liquid state 

A qualitatively new solution, however, is obtained 
within the present approximation for the disordered 
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FIG. 5: Kinetic energies E/t of the OL state for uncorrelated 
(U = 0, dashed line) and correlated (U = 00, full line) e g 
electrons (at 7 = 1), as functions of the electron density n. 



state, where double occupancies are on average elimi- 
nated by the slave bosons, and this correlation effect leads 
to an increase of the kinetic energy. The minimum en- 
ergy is obtained when the pseudofermion densities are 
equal, (fii+) = (ni-) = 5(1 — x), and the Gutzwiller 
renormalization factors take the simple form, 



= q±{x) 



2.1- 



(5.22) 



Then the pseudofermion bands, 

£g=co, ± (k) = «(s)ep=o 1± (k) 



(5.23) 



represent formally the superposition of the FO+ and 
AO± bands given by Eq. (|3.36j) . typical for uncorre- 
lated e g electrons, but now renormalized by correlations. 
They interpolate correctly between the case of uncorre- 
lated electrons in an empty band (x ~ 1) and a Mott 
insulator at half-filling (x = 0) where the dispersion is 
fully suppressed, as illustrated in Fig. [5] Owing to the 
Gutzwiller factors the kinetic energy has a minimum at 
filling n — 0.5, and approaches zero at n = 1. Thus, 
the kinetic energy has a similar doping dependence to 
that found in a spinless fermion model, i.e., for fermions 
with a single orbital flavor. As in the spin case^ one 
can argue that at x ~ strong correlations lead to an 
effective exclusion principle between the two degrees of 
freedom also in k space, i.e., for each momentum k only 
one orbital flavor may be occupied. 

This OL state is fully isotropic in the sense that the 
mean- field values of the pseudospin operators vanish, i.e., 



(If) = (If) = (If) = 0. 



(5.24) 



For the z-component this follows immediately from 
Eq. (|5.9[) once bf + = 6f_. For the other components 
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we apply Eqs. ([57X5]) and (f5T20)) to Eq. (|5TT0j) and obtain 



bi+bi- e 



-2ii), 



(fi-fi + ) 



(5.25) 



and similarly for (T^~ ). The pseudofermion averages 
can be determined by making use of Fourier transforma- 
tion: since the Fourier-transformed Hamiltonian (|5.21|) 
can be diagonalized analytically, the Fourier-transformed 
pseudofermion operators can be expressed in terms of the 
eigenvectors {e k +,e k _}, with the result 



ul + h-)+{fLU) = ^ 



(<4- 



e k+/-\ej, 



(/k+/k-) — (/k-/k+) — 



- '^v v k + k+/ 



(4- e k- 



(5.26) 

Since the eigenvalues £% l ^ 00 ± (k) are cubic invariant [see 
Eq. (|5.23[) ] in each of the two bands the three states with 
the components of k cyclically permuted are either all 
occupied or all unoccupied, and thus 



cos k. 



Si 



a i t \ 
-(e k+ e k+ ) 



(5.27) 



is independent of a, and similarly for C_. It then follows 
from the form of C k and L> k [see Eqs. (pH))) and ([XTU)) ] 

that the expressions (/ k ±/ k=F ), given by Eqs. (|5.26|) . both 
give zero when summed over the Brillouin zone, and so 



(ftf t -) = (flf i+ ) = 0, 



(5.28) 



and Eq. (jSTMl) follows. 

The absence of a preferred orientation of the pseu- 
dospin implies that there is no orbital preferentially oc- 
cupied. In particular, (Tf) = and (T?) — imply [see 
Eq. (E3)] that 



( C iz C iz 



J 

"ix ~-~tx I 



( C iz C ix d~ C ix C iz) ^, 



(5.29) 



from which it follows that the same relations hold for the 
operators {c' iz , c'\ x } obtained after an arbitrary U(l) ro- 
tation, as is easily verified explicitly or by observing that 
Tf and Tf rotate as an E doublet [compare the Ap- 
pendix]. Thus, the OL is SU(2) symmetric — random 
complex or random real orbitals are equivalent, and in- 
deed the identical OL state is obtained using real orbitals, 
as shown in the Appendix. This correlated disordered 
OL state with completely randomly occupied orbitals 
is apparently different from that proposed by Ishihara, 
Yamanaka, and Nagaosa^i in which the planar orbitals 
y 2 , y 2 — z 2 , z 2 — x 2 } play a prominent role. 



C. Absence of the Nagaoka theorem 



{x 2 



Before investigating the stability of the OL state in 
Sec. IV Dl let us consider the special case of a single 




FIG. 6: (Color online) Kinetic energies E per site at electron 
density n = 0.75 and U — oo for increasing off-diagonal hop- 
ping oc 7 in Eq. (|2.9p . as obtained in the KR approach for: 
the OL ground state (solid line), FO+ (dashed line), AO± 
state (AO, long-dashed line), and energy Eq (filled squares) 
for the ground state of a four-site plaquette (P). Also shown 
are the energies of the lowest two excited states for the pla- 
quette: a nondegenerate state which splits off the degenerate 
ground state at 7 = (empty squares), and a doubly degen- 
erate state with finite excitation energy at 7 = (diamonds) . 



hole in a half-filled system. In the spin case (7 = 0) the 
celebrated Nagaoka theorem^ one of the very few exact 
results in the theory of itinerant magnetism, then applies: 
Nagaoka has shown that the ground state is FM when 
a single hole/electron is added to a half- filled system, 
described by the spin Hubbard model at U = 00. A 
central assumption of this theorem is that the kinetic 
ener gy conserves the spin flavor (see, e.g., the proof in 
Ref. |48I ), precisely the feature not obeyed by the orbital 
flavor of e g electrons. Thus, at 7 7^ no exact statement 
can be made for the orbital Hubbard model (|2.9|) and, 
a priori, one expects that polarized states are harder to 
stabilize in this case. 

We have investigated the consequences of the SU(2) 
symmetry breaking, i.e. of the pseudospin non- 
conservation, by analyzing the exact solution for a pla- 
quette (four-site cluster) filled by three electrons, as a 
function of 7. In the spin model, at 7 = 0, the ground 
state, with kinetic energy E\j = — 0.25i per site, is four- 
fold degenerate, corresponding to maximum spin S = | 
as required by the Nagaoka theorem. At 7 > it splits 
into four nondegenerate states: the ground state and 
three excited states (the lowest of them is shown in Fig. 
[6]). The first excited state in the spin model (7 = 0) is 
doubly degenerate, and this degeneracy is not removed 
at 7 > 0, and the two states lower their energy when 7 
increases towards 7=1. For 7 > 0.4 this degenerate 
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excited state has already a lower energy than any other 
excited state (the level crossing is shown in Fig. [5]). None 
of these states can be classified by a pseudospin quantum 
number. In the genuine orbital case (7 = 1) the kinetic 
energy per site in the ground state, Eq — — 0.44i, is 
much lower than in the spin case (at 7 = 0), showing 
that a considerable amount of kinetic energy is gained 
when the orbitals get disordered and full advantage is 
taken of the pseudospin non-conserving hopping. This 
result suggests that a similar tendency towards disorder 
should be present in the thermodynamic limit. 



D. Stability of the orbital liquid phase 

Also for the full 3D model it is instructive to consider, 
at fixed density, the variation with 7 of the total energy 
E of possible ordered and disordered states. We do so 
in Fig. [6] at the same filling n = 0.75 as one has in the 
plaquette filled by three electrons, in order to enable a 
comparison with the exact results for that finite system. 
The energy of the polarized FO+ state does not depend 
on 7 [see Eq. (|3.31[) ]. while that of the AO± state follows 
from the dispersion given by Eq. (|3.32[) . and decreases 
linearly with 7. At 7 = 1 it comes very close to that of 
the FO+ state, but remains still a little bit higher. At 
7 = the polarized FO+ phase has a lower energy than 
the OL state, which confirms that FM states are stable 
in a range of filling close to n = 1 in the 3D Hubbard 
model42 The energy of the OL phase decreases gradu- 
ally with increasing 7, and becomes lower than that of 
the FO+ phase (which stays constant) at 7 ~ 0.25. It 
is remarkable that the energy decrease in the OL phase, 
when going from 7 = to 7 = 1, is quite large, and 
actually of similar magnitude as the exact result in the 
finite system. Hence, one finds that in spite of the renor- 
malization of the hopping by g(l/4) = 0.4, the (kinetic) 
energy in the OL state is substantially lower than in the 
AO± state. 

Next we consider the variation of the total energy E of 
ordered and disordered complex orbital states with elec- 
tron filling n (Fig. [7]). In the spin model (7 = 0) the FM 
phase has somewhat lower energy than the disordered 
OL state close to half- filling, in the range n > 2/3. 50 
Our approach reproduces in this limit the known result 
of the slave boson approach, which gives a FM ground 
state for any bipartite lattice with the density of states 
being an even function of energy4^ When 7 is increased, 
E-po does not change, whereas -Eao±i initially at zero for 
7 = 0, decreases cx 7, and at 7 = 1 surpasses the FO+ 
state at x = 0.27. Hence, the slave boson approach repro- 
duces here the result of the HF approximation for these 
states J£ However, in spite of the band narrowing oc q{x), 
which is appreciable at these electron densities near half- 
filling, considerably more (kinetic) energy is gained in 
the OL state. This is basically due to the fact that both 
hopping channels contribute, which gives rise to the large 
density of states over the full frequency range, and at 
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FIG. 7: (Color online) Kinetic energy E in the KR mean-field 
approximation as functions of n for: AO± (long-dashed lines) 
and OL state (full lines) for increasing 7 = 0, 0.5, 0.707, and 1 
from top to bottom; the dashed line shows the kinetic energy 
of the FO+ state which is independent of 7. 




Y 

FIG. 8: Region of stability of the FO states at U = 00 as a 
function of 7; the transition to the OL state from the FOx(z) 
and complex FO+ state are shown by the full and dashed line, 
respectively. 



small doping in particular [compare Fig.[TJd) with Fig. [2]. 
We may conclude that the presence of the additional non- 
pseudospin-conserving hopping channel, associated with 
the absence of SU(2) symmetry, implies that more kinetic 
energy can be gained by paying correlation energy than 
in the spin case, and that this favors the disordered OL 
state sufficiently to make its energy lower than those of 
the complex orbital-ordered states at any value of n. 
Finally we compare at U = 00 the energies of all states, 
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both with complex and real orbitals, varying n and 7. 
One finds that AO states are never stable in this limit 
of strong correlation, while FO states are stable only at 
small 7 [Fig. [S]. At 7 = (the spin case) the FO+ and 
FOx (FOz) states are necessarily degenerate, but at any 
7 > the phases with ordered real orbitals have lower 
energy, with FOz (FOir) being more stable at n < 0.71 
(n > 0.71). The range of FO order shrinks gradually 
with increasing 7, and above 7 ~ 0.94 the OL phase is 
stable in the entire range of n. We argue that at finite 
U the kinetic energy will become even more dominant 
and thus will strongly favor disorder, except near n ~ 
1 where SE stabilizes real-orbital AO order i2iMiSi2£ We 
thus conclude that for the e g orbital Hubbard model (7 = 
1) doping triggers a crossover to the OL state at any 
U , supporting earlier conjectures that such a disordered 
state is realized . 47 i 51 



E. Brinkman-Rice transition at n — 1 

At half- filling (n = 1) it is straightforward to apply the 
finite-?/ version of the KR formalism^ and investigate 
the generic metal-insulator transition in the orbital dis- 
ordered phase, ignoring the AO order promoted by the 
SE. Here one introduces as a counterpart to the bosons 
e, which control the empty configurations |0) = ej|vac), 
also bosons di which control the double occupancies 
<4|-cjjj0) = d\f^fj, I vac). The mean- field approximation 
gives then the renormalization factor (at n = l), 14 

77(d) = 8d 2 (l - 2d 2 ), (5.30) 

where d = (di) is the average amplitude of a doubly 
occupied configuration in the ground state. The bands 
are then given by the dispersion for free electrons (|3.35|) 
renormalized by 77(d), 

e%% (k) - v{d)e u= o,± (k) = -r]{d)t[A k ± 7 5 k ] . (5.31) 

So the kinetic energy is 77(^)60(7), where 60(7) is the ki- 
netic energy of the uncorrelated OL, obtained by inte- 
grating the two bands £a=o,±(k) (|3.35[) up to half-filling, 
while the Coulomb repulsion gives an energy Ud 2 per 
site. 

For the spin model (7 = 0) this problem was solved 
by Brinkman and Rice^ who showed that an 'insulat- 
ing' state (with d = 0) sets in above U c ~ 8t (in the 
present units). It is well understood by now (see Ref. |22[) 
that this mean-field theory does not give an accurate de- 
scription of the metal-insulator transition (in particular 
it ignores all charge fluctuations in the insulating phase, 
where in reality d ^ 0)^ By analogy, one expects that 
also in the present case d ^ at any finite U, and in fact 
this follows from the large-?/ expansion analyzed for the 
ordered phases in Sees. IIV Bl and IIV CI (for a disordered 
phase a similar analysis could also be made). Neverthe- 
less, the Brinkman-Rice transition from a 'metallic' to 
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FIG. 9: Metal (M) to insulator (I) transition with increasing 
U, as obtained in the disordered phase as a function of 7 at 
n — 1. At 7 = 1 (orbital case) the transition occurs in the 
ground state, while at 7 = (spin case) the result of Ref. [13 
is reproduced^ 

an 'insulating' state at n — 1 illustrates nicely the com- 
petition between kinetic energy and Coulomb repulsion 
energy^ 4 - and so it is worthwhile to consider the gen- 
eral case, i.e., with arbitrary 7. Then, completely anal- 
ogously to the spin case, an 'insulating' state is found 
above U c (j) — 8|eo (t) I - Similar to what happens upon 
doping (i.e., at finite x) in the U = 00 limit consid- 
ered above, here upon allowing double occupancy (i.e., 
finite d) at n = 1, the metallic phase gains additional ki- 
netic energy oc 7 due to non-pseudospin-conserving hop- 
ping which lowers the kinetic energy below the value due 
to pseudospin-conserving hopping alone (the only one 
present in the spin case). Therefore, the metallic phase 
survives up to a higher value of U than at 7 = 0, as 
shown in Fig. [5] 

VI. SUMMARY AND CONCLUSIONS 

In this paper we have made a detailed analysis of the 
e g -orbital Hubbard model on a cubic lattice, exploring 
the consequences of the absence of SU(2) symmetry and 
highlighting them by making a comparison with the fa- 
miliar SU(2)-symmetric spin Hubbard model. In the first 
part we studied the orbital-ordered phases, of which there 
is a great variety, precisely because of the lower sym- 
metry, emphasizing the difference between the complex- 
orbital states which retain cubic symmetry, and the real- 
orbital states in which cubic symmetry is broken. Ana- 
lytical results for the order parameter and the energy of 
each of these phases in the HF approximation at large 
U/t were presented, demonstrating that the total energy 
can be conveniently divided into two contributions: a ki- 
netic energy oc t given by the U = 00 limit, and a SE 
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contribution oc t 2 /U . The SE decides about the relative 
stability of the various phases at half-filling, while the 
kinetic energy contributes and finally becomes dominant 
upon doping. This analytical treatment allowed us: (i) 
to demonstrate explicitly that SE contributes in both AO 
and FO states, (ii) to demonstrate that the real-orbital 
states have their orbitals aligned with the cubic axes, as 
well as (iii) to elucidate the structure of the HF phase dia- 
gram for the ordered phases obtained numerically ! 20 ' 23 ' 24 
We emphasize that these properties of e g orbital degrees 
of freedom are essentially different from those of ti g ones, 
because the latter satisfy certain symmetries and are thus 
conserved in the hopping processes^i 5 - 5 - 

In the second part we investigated the disordered 
orbital-liquid state. We have demonstrated that in the 
strong-correlation limit (U ^S> t) indeed orbital (FO or 
AO) order is not robust for e g orbitals, and gets replaced 
by a disordered (OL) phase, if one goes beyond the HF 
approximation and includes electron correlation effects in 
the disordered phase as well. This leads us to the con- 
clusion that the HF results ; 20 ! 23 ' 24 suggesting that either 
the FO+ or the AO± state is realized in a broad range of 
doping, are particularly misleading for the orbital Hub- 
bard model. Here the present findings agree qualitatively 
with the results of the self-consistent second order pertur- 
bation theory obtained by Kubo and Hirashima^ The 
situation could be somewhat different in the 2D case, 
however, where a tendency towards particular orbital or- 
derings with larger amplitude of x 2 — y 2 orbitals is favored 
by geometry. 28 ' 38 

We considered specifically the U = oo limit, where the 
OL competes with fully polarized ordered phases and we 
have shown that it is more stable than any of either uni- 
form FO (|4.21|) or staggered AO (|4.22j) states. However, 
at finite U and for sufficiently low doping x, real-orbital 
C- AO order is stabilized by a superposition of the SE and 
the JT effect. Particularly in the regime of low doping 
the JT interactions might be stronger than the electronic 
interactions of double-exchange type, and the induced or- 
bital order dictates then the type of magnetic orde n 35 ' 57 
This regime is particularly difficult in realistic models for 
manganites, as the orbital interactions induced by oxygen 
distortions^ and the orbital polarization around doped 
holes 58 give additional important contributions and sup- 
port particular types of orbital order. Furthermore, the 
overall stability of ordered versus disordered (OL) phases 
changes when a realistic Hund's coupling is included. 59 
It has been shown that the FM phase shrinks then to a 
range of doping 0.2 < x < 0.5, the A-type AF phase is 
stable near x = 0.5, while the C-AF phase takes over at 
higher hole doping. 

Summarizing, the absence of SU(2) symmetry in the 
e 9 -orbital Hubbard model has severe consequences for 
the properties of the model itself and for the stability of 
orbital-ordered states. The Nagaoka theorem does not 
apply to the model of correlated e g electrons at U = oo, 
ordered states are harder to realize than in the spin 
case, and the Brinkman-Rice transition occurs at a higher 



value of U. The qualitatively different properties of the 
ordered phases show up most clearly in the inverted sta- 
bility (with respect to the spin case) of the ordered phases 
with complex orbitals, with ferro (staggered) orbital or- 
der favored at small (large) doping. Most importantly, 
the exciting suggestion that such complex-orbital ordered 
states could be stable at finite dopin g 20 ' 23 ' 24 has been 
disproved, because of the inherent tendency of e g sys- 
tems towards orbital disorder due to the enhancement of 
the kinetic energy when SU(2) symmetry is absent. All 
these features show that several properties of spin sys- 
tems which are usually taken for granted, such as: (i) 
the very fact that a ferromagnetic state is an eigenstate 
of either an itinerant or the Heisenberg Hamiltonian, and 
(ii) the absence of superexchange in ferromagnetic states 
— are in fact the consequences of the SU(2) symmetry 
of the respective spin models. 
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APPENDIX: SLAVE BOSON 
REPRESENTATION FOR REAL ORBITALS 

The real-orbital version of the transformation of the 
electron operators to slave boson and pseudofermion op- 
erators may be derived by making repeated use of the 
relation between the real and complex orbitals, as given 
by Eqs. (|2.4[) . Thus with the real-orbital electron opera- 
tors given by 

4 = 75 (4 + 4-): 

4 = ^(c\ + -c\_), (A.l) 

we similarly define the real-orbital pseudofermion opera- 
tors by 

ft = 75 (/*+ + /<-)> 
ft = js(fl-fl), (A.2) 
while for the slave boson operators we set 

4 = ^(4 + 4), 

4 = 7l(>l -4)- ( A - 3 ) 

Then the fermions (electrons) transform under U(l) ro- 
tations as 

Ui(e)4 z Uj(e) = cos(tf/2) C L-sin(0/2)4, 
^(0)4^(0) = Hk(0/2)4 + cos(0/2)4, (A.4) 
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and similarly for the pseudofermions, while the slave 
bosons transform as 



is readily verified to eliminate the unphysical singly- 
occupied states, 



U t (9)blul(9) = - aindbl + co S 0b\ x . (A.5) 



sin 9 b. 



t 

ix ' 

t 



The different sign choice in Eq. (|A.3j) as compared to 
Eqs. (jA.lj) and ()A.2j) makes the slave bosons rotate in the 
opposite direction as the (pseudo)fermions. This com- 
pensates for the doubled rotation angle in the sense that 
the transformations are identical for slave bosons and 
(pseudo)fermions when 9 is a cubic angle, and so the 
pairs {c\ z ,c\ x }, {fj z ,fj x }, and {b\ z ,b\ x } all transform as 
the 9 and e component of a cubic E doublet. 

Substituting the complex-orbital slave boson represen- 
tation (|5.4[) into Eq. (|A.1[) and applying the inverse trans- 
formations to (|A.2|) and (|A.3[) . one obtains the slave bo- 
son representation for the real-orbital fermionic operators 
{ c L' c Li analogous to Eq. (|5.4[) . The result is 

c 1 - +^-lV f 1 - 6 f f 1 V 

4 = -^(blfL+blfl)^, (A.6) 

corresponding to a representation of the local states by 

| ? 0) = et|vac), 
I«>=4li0) = +^(4.4-4/l)|vac), 

-^ 2 {b\jl+b\ z fl)\^c). (A.7) 



iO) = -- 



One recognizes that Eqs. (|A.6p are indeed the proper 
expressions for the E doublet resulting from the prod- 
uct representation E (g) E (g) A\ ^ The expressions (|A.6|) 
are actually even U(l)-invariant, i.e., after a rota- 
tion in orbital space by an arbitrary angle 9, they 
also hold between the fermion operators {c%,c'l,} — 

{Ui(8)4 z Uj (8) ,Ui(0)4 x Uj (9)}, transformed according 
to Eq. (|A.4j) , and the slave boson and pseudofermion op- 
erators {b'i Z , b'^} and {f' iz , f' ix }, transformed according to 
Eqs. (|A.5|) and (|A.4[) . respectively. Consequently, since 
the hopping Hamiltonian (|2.3|) is invariant under a trans- 
formation (|A.4jl of the fermion (electron) operators when 
9 is one of the cubic angles 0, ±47r/3 and is accompanied 
by the corresponding permutation of the cubic axes, this 
cubic invariance is retained when the Hamiltonian is ex- 
pressed in terms of the slave boson and pseudofermion 
operators by means of Eq. (|A.6|) . 

The constraints given by Eqs. (|5.7[) are now replaced 

by 



b \Az + b\ x b lx 
b\Ax ~ h \Az 



4 e i 



ft f. i ft f. 



fl f. - f! f 

J IZJ IX J IX J 'I 



(A. 



Again the first constraint excludes double-occupancy, as 
required in the limit U = oo, while the last constraint 



= ^ 2 {b\jl + b\J 



vac) 



IZ •> IX 



(A.9) 



When the constraints are obeyed rigorously and the un- 
physical states strictly projected out, operators connect- 
ing the physical and unphysical subspaccs necessarily 
vanish identically. Specifically one finds 



fct b- - 6 f b- 

^%Z IZ ^ IX IX 



ft f. 

J izJ IZ 



ft f. 

J ix J i: 



o, 



b\zKx + b\ x b lz = fjj lx + fjj lz = 0. (A.10) 

It is obvious from the above that the earlier attempt 
made in Ref. [fill to construct a real-orbital slave bo- 
son representation by means of c\ z — b\ z j\ z ei and c\ x — 
b\ /t 

IX J IX 

factors by 



b\ x f ix ei 1 followed by renormalization of the slave boson 



vt 



6 f e- 



\f(l-4e t -bib lx )(l-bib t S 



6 f e- 



V /(l-ete i -6t z& .J(l-6t ;c6 .j 



, (A.ll) 



was misguided because it does not conserve the cubic 
symmetry, and is thus bound to lead to spurious re- 
sults. However, also the present real-orbital representa- 
tion, though invariant in itself, leaves us with the prob- 
lem to construct a proper cubic-invariant renormaliza- 
tion. This is not straighforward because the hopping 
Hamiltonian, when expressed completely in terms of slave 
boson and pseudofermion operators referring to l z' and 
l x\ takes a different appearance for each cubic axis, like 
in Eq. (|2.3p . Moreover, the apparently plausible renor- 
malization by means of Eqs. (|A. 1 1|) is not allowed even 
in combination with the representation (|A.6|) . because z\ z 
and z\ x as defined by Eqs. (|A. 1 1|) do not constitute a cu- 
bic E doublet as their denominators are not cubic invari- 
ants. Having them replace b\ z e i and b ix e i in Eqs. (|A.6p 
would spoil also the cubic E doublet nature of the thus 
renormalized c\ z and cj x , and so destroy the cubic sym- 
metry of the Hamiltonian. Equally seriously, it would 
also cause the Hamiltonian to commute no longer with 
the constraints. 

A renormalization not suffering from the above prob- 
lems and still in the spirit of the Kotliar-Ruckenstein 
AnsatzW is given by 



U lZ C l 



(i-ek-§nf)(i-§nf>) 



b\ x e-i 



(A.12) 
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where n\ = b\ z b iz + b\ x b ix . The mean- field approxima- 
tion is now made, as in Sec. IV Al by replacing only the 
amplitudes but not the phases by c-numbers. So, for the 
offdiagonal two-boson products we set, similarly to what 
was done in Eqs. (|5.15|) . 

( b L e i) = (4 b iz) = hei cos(a t - 

{b\ x e t ) = (e\b ix ) = b t e % sinfo, - ^). (A.13) 

where bi and Si are again real quantities. For the diagonal 
two-boson products we set 

Me,) = el (A.14) 

Actually, the real-orbital boson occupation numbers 
n\g = b\ z b iz and nf x = b\ x b ix are not invariants with 
respect to U(l) rotations, and so one would prefer to 
set, in accordance with Eqs. (|A.13j) . the corresponding 
diagonal averages equal to 

(b\ z b lz ) = 11 cos'te - 

(b\ x b lx ) = b\ S in 2 (a t -A), (A.15) 

in order to make them transform in the same way as the 
occupation numbers, by setting also 

(b\A x ) = (b\A z ) ee \b\ sin(2a, - 2^). (A.16) 

However, the expressions (jA~T5|) and (fA~T6|l do not sat- 
isfy Eqs. (|A.10p . and so it appears to be impossible to 
assign a nontrivial dependence on the phase operator i?i 
to {b\ z b iz ) and (b] x b ix ) and yet simultaneously respect 

(Oof . 



The issue is immaterial for carrying out the KR pro- 
cedure, since in a state with uniform density, i.e. with 
ef = x for all i, it follows from the constraints (|A.8j) that 
both for (|A.14[) and for (|A.15[) the amplitude satisfies 
bf = 1 — x, so that 

(4) = (hz) = VMx) cos(a, - Si), 

(4) = {z ix ) = ^/Hx) sm{ ai -k). (A.17) 

Inserting this into Eqs. (|A.6|) and defining new pseud- 
ofermions by 

ft = cos a t fj z + sin a, f\ x , 

= cos(i?i - en) f} z +sin(i?i - oti) fj x , 

ft = ~ sin a = ft + cos °H ft i 

= sin(i?i - a*) fj z - cos^j - on) ft, 

(A.18) 

where {ft, ft} are related to {// + , f}_} [sec Eq. PJ] 
by Eqs. (|A.2[) . one finds that the mean-field approxima- 
tion effectively leads to the replacements 

4 = VoWfl, 4 = V^/L . (A.19) 

The kinetic part of the Hamiltonian is thus simply renor- 
malized by the Gutzwiller factor q(x), exactly the same 
result as obtained in the complex-orbital approach. As 
the Hamiltonian is therefore again cubic, it follows that 
the resulting real-orbital OL is isotropic (i.e. ai = 7r/4 at 
all sites, and (b\ z b lz ) = (b\ x b ix ) = (1 - x)/2), and identi- 
cal to the OL obtained in the complex-orbital approach. 
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